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A semiclassical model based on the solution of the Vlasov equation for finite sys- 
tems with a sharp moving surface has been used to study the isoscalar quadrupole 
and octupole collective modes in heavy spherical nuclei. Within this model, a 
unified description of both low-energy surface modes and higher-energy giant res- 
onances has been achieved by introducing a coupling between surface vibrations 
and the motion of single nucleons. Analytical expressions for the collective re- 
sponse functions of different multipolarity can be derived by using a separable 
approximation for the residual interaction between nucleons. The response func- 
tions obtained in this way give a good qualitative description of the quadrupole 
and octupole response in heavy nuclei. Although shell effects are not explicitly 
included in the theory, our semiclassical response functions are very similar to the 
quantum ones. This happens because of the well known close relation between 
classical trajectories and shell structure. The role played by particular nucleon 
trajectories and their connection with various features of the nuclear response is 
displayed most clearly in the present approach, we discuss in some detail the damp- 
ing of low-energy octupole vibrations and give an explicit expression showing that 
only nucleons moving on triangular orbits can contribute to this damping. 



1. Introduction 

It is well known that the isoscalar quadrupole and octupole response of nu- 
clei displays both low- and high-energy collective modes ^. Also known is 
that semiclassical models have difficulties in describing both these system- 
atic features of the isoscalar response, in particular, models based on fluid 
dynamics, see e.g. ^, can explain the giant resonances, but fail to describe 
the low-energy collective modes. On the other hand it is known from quan- 
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turn studies that the coupling between the motion of individual nucleons 
and surface vibrations plays an essential role in low-energy nuclear collec- 
tive modes, see e.g. Semiclassical models of the fluid-dynamical type 
do not contain explicitly the single-particle degrees of freedom, so they can 
not describe the coupling between individual nucleons and surface motion. 

In the present contribution we review a study the isoscalar collective 
modes in nuclei ^'"^ made by using a semiclassical approach that includes 
the single-particle degrees of freedom explicitly and thus allows for an ac- 
count the coupling between individual nucleons and surface motion. Our 
model is based on the linearized Vlasov kinetic equation for finite systems 
with moving surface The coupling between the motion of individual 
nucleons and the surface vibrations is obtained by treating the nuclear sur- 
face as a collective dynamical variable, like in the liquid drop model. Here 
we concentrate our attention on the isoscalar quadrupole and octupole col- 
lective modes in heavy spherical nuclei, an application of the same model 
to the compression dipole modes has been discussed in a previous meeting 
of this series 

2. Reminder of formalism 

This Section recalls briefly the formalism of References ^'^ which is at the 
basis of the present approach. 

The fluctuations of the phase-space density induced by a weak external 
force can be described by the linearized Vlasov equation, which is usually a 
differential equation in seven variables. For spherical systems this equation 
can be reduced to a system of two (coupled) differential equations in the 
radial coordinate alone ^ . This is achieved by means of a change of variables 
and a partial-wave expansion: 



The functions 5fj^jq{e, A, r, w) are partial- wave components of the (Fourier 
transformed in time) density fluctuations for particles with energy e, magni- 
tude of angular momentum A and radial position r, the ± sign distinguishes 
between particles having positive or negative components of the radial mo- 
mentum Pr . 

The other terms in the expansion are Wigner matrices and spherical 
harmonics. 



5f{v, ip,uj)= [S/m'n (e. A, r, lu) + Sf^-^ (e. A, r, lu)] 



LAIN 




(1) 
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In order to solve the one-dimensional linearized Vlasov equation for the 
^fMN functions we must specify the boundary conditions satisfied by these 
functions. Different boundary conditions allow us to study different physical 
properties of the system, so the fixed-surface boundary conditions employed 
in ^ were adequate to study giant resonances, but different (moving-surface) 
boundary conditions ^ must be introduced in order to study surface modes. 
We assume that the external force can also induce oscillations of the system 
surface according to the usual liquid-drop model expression 

R{i», ^,t) = R + Y^ SRLM{t)YLM{^, (2) 

LM 

and the boundary condition satisfied by the functions 5f^^ at the nuclear 
surface is taken as 

5fMN{R) - 5fMN{R) = 2F\e)lUjpr5RLM{u^) . (3) 

This equation has been derived with the assumption that the equilibrium 
phase-space density is a function i^(e) of the particle energy alone, F'(e) 
is its derivative. The boundary condition (3) corresponds to a mirror re- 
flection of particles in the reference frame of the moving nuclear surface, it 
provides a coupling between the motion of nucleons and the surface vibra- 
tions. A self-consistency condition involving the nuclear surface tension is 
then used to determine the time (or frequency) dependence of the additional 
collective variables 5Rlm{^) ■ 

Now, assuming a simplified residual interaction of separable form, 

v{ri,r2) = ULr^r^ , (4) 

the moving-surface isoscalar collective response function of a spherical nu- 
cleus, described as a system of A interacting nucleons contained in a cavity 
of equilibrium radius R = 1.2^3 fm, is given by 

nL{s)^nL{s) + SL{s). (5) 

Instead of the frequency w, as independent variable we have used the more 
convenient dimensionless quantity s = uj/ {vp/R) {vf is the Fermi velocity). 
The response function TZl{s), given by 

describes the collective response in the fixed-surface limit. The response 
function 72.^ (s) is analogous to the quantum single-particle response func- 
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tion and it is given explicitly by 

Kis) = E E (^^^)' /' 'r.N{x)-^f^^^^ , (7) 

where is the Fermi energy and the quantity £ is a vanishingly small 
parameter that determines the integration path at poles. 
The functions s„jv(a;) are defined as 

, , riTT + arcsinfx) 

X 

The variable x is related to the classical nucleon angular momentum A. 

The quantities Cln in Eq. (7) arc classical limits of the Clcbsh-Gordan 
coefficients coming from the angular integration. In principle the integer N 
takes values between —L and L, however only the coefficients Cljv where 
N has the same parity as L arc nonvanishing. The coefficients qI^^{x) 
appearing in the numerator of Eq. (7) have been defined in Rcf. ^, they are 
essentially the classical limit of the radial matrix elements of the multipole 
operator and can be evaluted analytically for L = 2, 3. 

The function Sl{s) in Eq. (5) gives the moving-surface contribution to 
the response. With the simple interaction (4) this function can be evaluated 
explicitly as ^'^ 



Sl{s) = -. 



1 - KL-R-lis) [Cl - Xl(s)][1 - klTZHs)] + klR^xUs) + QoR^^V 

(9) 

with Cl = aR^{L - f)(L + 2) + {Cl)couI, o « IMeVfm-^ is the sur- 
face tension parameter obtained from the mass formula, {Cl)couI gives the 
Coulomb contribution to the restoring force and go = A/^R^ is the equi- 
librium density. The functions Xl{^) ^^'^ Xl{s) are given by 



Xl(s) = E E ((^LNf I dxx^ s„jv(a;)- 



n=-ooiV=-L S + ie-SnN{x) 

(10) 

and 



Xl{s) = -—eF{s + ie)y^{CLNf' I dxx 



S + ie - SnN^X) ' 



(11) 



their structure is similar to that of the zero-order propagator (7). 

We refer to the papers ^'^ for further details on the formalism and discuss 
here only the main points. 
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Equation (9) is the main result in the present context. Together with 
Eqs.(5) and (6), this equation gives a unified expression of the isoscalar 
response function, including both the low- and high-energy collective exci- 
tations. By comparing the fixed- and moving-surface response functions, 
we can appreciate the effects due to the coupling between the motion of 
individual nucleons and the surface vibrations. 



3. Fixed- vs. moving-surface strength distributions 

The strength function Sl{E) associated with the response function (5) is 
defined as {E = hco) 

Sl{E) = --ImTlLiE) . (12) 

TT 

We discuss here the isoscalar quadrupole and octupolc strength distribu- 
tions. The strength Kl of the residual interaction (4) can be estimated in 
a self-consistent way, giving (^^, p. 557), 

with the parameter uiq given by ojq « 41A~i'McV . Since this estimate is 
based on a harmonic oscillator mean field and we are assuming a square- 
well potential instead, we expect some differences. Hence we determine the 
parameter kl phenomenologically, by requiring that the peak of the high- 
energy resonance agrees with the experimental value of the giant multipole 
resonance energy. This requirement implies kl ~ 2kbm for L = 2,3. 

In Fig.l we display the quadrupole strength function (L=2 in Eq. (12)) 
obtained for A = 208 nucleons by using different approximations. The 
dotted curve is obtained from the zero-order response function (7), it is 
similar to the quantum response evaluated in the Hartree-Fock approx- 
imation. The dashed curve is obtained from the collective fixcid-surface 
response function (6). Comparison with the dotted curve clearly shows the 
effects of collectivity. The collective fixed-surface response has one giant 
quadrupole peak. Our result for this peak is very similar to that of the 
recent random-pasc approximation (RPA) calculations of (cf. Fig. 5 of 
^^). However, contrary to the RPA calculations, there is no signal of a low- 
energy peak in the fixed-surface response function. The solid curve instead 
shows the moving-surface response given by Eqs. (5) and (9). Now a broad 
bump appears in the low-energy part of the response and a narrower peak 
is still present at the giant resonance energy. Of course the details of the 
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low-energy excitations are determined by quantum effects, nonetheless the 
present semiclassical approach does reproduce the average behaviour of this 
systematic feature of the quadrupole response. 

We finally notice that the width of the giant quadrupole resonance is 
underestimated by our approach, this is a well known limit of all mean-field 
calculations that include only Landau damping. A more realistic estimate 
of the giant-resonance width would require including a collision term into 
our kinetic equation. 
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Figure 1. Quadrupole strength function for a liypotiietical nucleus of A = 208 nucleons. 
The dotted curve shows the zero-order aproximation, the dashed curve instead shows the 
collective response evaluated in the fixed-surface approximation. The full curve gives the 
moving-surface response. 



In Fig. 2 we show the octupole strength function (L = 3 in Eq. (12)). 
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The zero-order octupole strength function (dotted curve) is concentrated 
in two regions around 8 and 24 McV. In this respect our semiclassical re- 
sponse is strikingly similar to the quantum response, which is concentrated 
in the Ifuv and 3tiu) regions. This concentration of strength is quite re- 
markable because our equilibrium phase-space density, which is taken to be 
of the Thomas-Fermi type, does not include any shell effect, however we 
still obtain a strength distribution that is very similar to the one usually 
interpreted in terms of transitions between different shells. 

We can clearly see that the collective fixed-surface response given by 
Eq. (6) (dashed curve) has two sharp peaks around 20 Mev and 6-7 Mev. 
The experimentally observed ^ concentration of isoscalar octupole strength 
in the two regions usually denoted by HEOR (high energy octupole reso- 
nance) and LEOR (low energy octupole resonance) is qualitatively repro- 
duced, however the considerable strength experimentally observed at lower 
energy (low-lying collective states) is absent from the fixed-surface response 
function. The most relevant change induced by the moving surface (solid 
curve in Fig. 2) is the large double hump appearing at low energy. This 
feature is in qualitative agreement both with experiment ^ and with the 
result of RPA-type calculations (see e.g. ^^). We interpret this low-energy 
double hump as a superposition of surface vibration and LEOR. 

The moving-surface octupole response of Fig. 2 displays also a novel 
resonance-like structure between the LEOR and the HEOR (at about 13 
MeV for a system oi A = 208 nucleons) . 

We also find that the parameters dRsMit), describing the octupole 
surface vibrations in Eq. (2), approximately satisfy an equation of motion 
of the damped oscillator kind: 

Ds5RsM{t) + 73SR3M{t) + C^SRsMit) = . (14) 

The friction coefficient can be evaluated analyticallty in the low- 
frequency limit, giving (for a generic L) 

IL = Iwf ^ 2L + 1 ^ N ' ^^^^2 ' 2^ ' 2^cosa„jvsm a„jve(--Q„jv) , 

N=l n=l 

(15) 

with 7^/ = jQoppR'^ and a„Ar = ^tt. The angles a„jv are related to 
the nucleon trajectories. In the octupole case the coefScient 7l=3 gets a 
contribution only from the term with n = 1 and N = 3, thus we see that 
only nucleons moving along closed triangular trajectories can contribute to 
the damping of surface octupole vibrations. 
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Figure 2. The same as in Fig.l for the octupoie strength function. 



4. Conclusions 

A unified description of the low- and high-energy isoscalar collective 
quadrupole and octupoie response has been achieved by using appropri- 
ate boundary conditions for the fluctuations of the phase-space density de- 
scribed by the linearized Vlasov equation. The response functions obtained 
in this way give a good qualitative description of all the main features of the 
isoscalar response in heavy nuclei, i. e. low-lying quadrupole and octupoie 
collective modes, plus quadrupole and octupoie giant resonances. In our 
model the low-energy modes arc surface oscillations and the coupling be- 
tween single-particle motion and surface vibrations is described by simple 
analytical expressions. 
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